"0","```r
?estimateR # read more about this function, and how ACE and Chao1 are calculated
alt_richness <- estimateR(scpdSM)
# Bind all richness estimates together and compare:
richness_all <- rbind(alt_richness, richness)
# plot them against one another
y.max <- floor(max(richness_all))+1 # set y max for plotting 
mainEstimates <- c('S.obs', 'S.chao1', 'S.ACE', 'S') 
plotRows <- match(mainEstimates, rownames(richness_all))
plotCols <- c(\"black\", \"blue\", \"green\", \"red\")
  
i=1
plot(richness_all[plotRows[i],] ~ ages, type=\"l\", ylim=c(0, y.max), xlim=c(21000, 0), xlab=\"Age (cal BP)\", ylab=\"Richness\", col=plotCols[i])
for (i in 2:length(plotRows)){
    lines(richness_all[plotRows[i],] ~ ages, col=plotCols[i])
}
```"
"0","```r
legend(\"bottomleft\", legend=mainEstimates, lty=1, col=plotCols, bty=\"n\", cex=0.75)
```"
